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(54) Process and apparatus for controlling the growth of a crystal 

(57) The invention concerns a process and an 
apparatus for controlling the growth of a crystal, which 
growth is governed by a set of measurable and non- 
measurable variables. The process comprises the fol- 
lowing steps: 

establishing an on-line simulation software working 
with a reduced number of variables, the reduction of 
variables being performed by using a projection 
algorithm: 

speeding up the on-line simulation software by gen- 
erating data banks in which values of off-line pre- 
calculated variables are stored; 
tuning the on-line simulation software by adjusting 
the results predicted by on-line simulations to the 
results obtained by off-line simulations and by 
measurements; and 

establishing a control loop and controlling at least 
one of the variables in real time, the control loop 
using the speeded up and tuned on-line simulation 
software as an on-line observer. 
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Description 

Technical Field 

5 The invention is related to the design of a complete multivariable predictive control system for the regulation of sin- 

gle aystal growth processes. 

Backoround Of The Invention 

10 The growth of large single crystals of a specified quality is an important issue in many high-tech applications, such 
as electronic devices, fiber optics converters, lasers or infrared lenses. While many processes and variants can be used 
for this purpose, the most common way to grow large crystal boules is to resort to the Czochralski, Roat Zone, Horizon- 
tal Bridgman, or Vertical Bridgman technique. With regard to the production of crystalline films, epitaxy is the most com- 
mon technique. 

15 

Czochralski growth: 

In the first case, which represents the most frequently used apparatus, a heater (often an ohmic resistor) radiates 
power towards a crucible containing the molten material, hence maintained above the solidification temperature. 

20 A liner may separate the crucible from the melt. In the case of liquid encapsulated growth, a particular encapsulant cov- 
ers the melt to prevent evaporation of volatile components. An alternative method to supply energy to the crucible and 
the melt (which is frequently used to grow oxide crystals) is to resort to an induction heating system, in which applying 
a strong alternating current to the coil generates heat directly within the crucible and/or the melt. 

At the outset of the process, a crystal seed is dipped into the liquid and then slowly pulled while a cylindrical crystal 

25 is progressively grown. The crystal radiates heat towards the surrounding environment, which consists of the outer shell 
of the furnace (cooled by water circulation) and. in some cases, of a high pressurized gas. In this latter situation, the 
release of heat by conduction and convection in the gas may be important. It is in practice necessary to let both crystal 
and crucible rotate around the vertical axis, in order to ensure axisymmetry of the overall process, and to control the 
buoyancy-driven flow in the liquid. 

30 The solidification front between crystal and melt is the place where interactions between the liquid and solid phases 
take place. A major difficulty of growing high quality crystals often lies in the dependence of interface shape and segre- 
gation upon the melt flow pattern. Melt convection can be induced by temperature or solute gradients within the liquid 
(natural or thermosolutal convection), by crystal and crucible rotations (forced convection), or by surface tension gradi- 
ent in the liquid-environment interface (Marangoni convection). In semiconductor growth, the Grashof number of the 

35 flow is fairly high, which may give rise to an oscillating or even turbulent behavior in the melt. These effects may have 
severe drawbacks, such as crystal striations, and a solution to reduce their importance is to generate a strong magnetic 
field, which plays a role as a vigorous brake in the flow. 

The growth of oxide crystals presents some important differences with respect to semiconductors. While typical 
Grashof numbers are lower, heat transfer in the melt becomes convection dominated in view of the higher Prandtl 

40 number of molten oxides. Another peculiar feature is the importance of internal radiation, which may be a significant 
mode of heat transfer in the crystal. 

Bridgman growth: 

45 An alternative solution is to resort to the horizontal or vertical Bridgman technique. Here, the molten material is con- 
tained in a cylindrical ampoule, at the bottom of which solidification starts from a crystal seed (after back-melting). The 
outer thermal environment of the ampoule consists of a hot. an insulation and a cold zone. During growth, the ampoule 
is slowly shifted downwards, and the temperature distribution induced by the outer ertvironment in the solid and molten 
material shifts upwards accordingly, white keeping more or less the same axial profile. Hence, melt solidifies progres- 

50 sively until the crystal is completely formed. 

A major advantage of vertical Bridgman growth is that the solid phase is located under the liquid phase and, con- 
sequently, that the temperature gradient is pointing essentially upwards in the melt. If isotherms were exactly horizontal, 
the layered fluid would be perfectly stable with respect to thermal convection. The problems related to the onset of a 
buoyancy-driven flow in the melt are, therefore, strongly reduced (a destabilizing radial temperature gradient is however 

55 generally present). The situation is more complex when considering solutal convection, which may be important for a 
non-dilute alloy Indeed, there is no general rule concerning the dependence of density with respect to solute concen- 
tration, and there may be solute rejection or incorporation at the interface. Hence, the solutal field can have a priori a 
stabilizing or destabilizing effect. It is worth noting that melt convection outside the diffusion-controlled boundary layer 
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can have a beneficial mixing influence on crystal uniformity. 

An efficient way to control the temperature gradient at the outer surface of the crystal and the melt is to replace the 
classical Bridgman apparatus by a dynamic gradient freeze apparatus. Basically, the heating system consists of a set 
of controlled annular elements monitored by appropriate thermocouples, which are placed mostly around the periphery 
5 of the growth vessel. It is no longer necessary to let the ampoule move downwards in the vessel, since the evolution of 
the outer temperature profile Is achieved by a slow vertical shift of power supply in the heating elements. 

Floating zone growth: 

10 At last, in the Float Zone technique, the material is in the form of a free-standing rod clamped only at its ends, in 
which a small zone is melted by suitable heating equipment. The melt is suspended like a drop between the two parts 
of the rod. The molten zone is moved through the rod over its whole length by heater or rod motion. A single crystal can 
be generated by spontaneous nucleation. or using a single-crystalline seed crystal as the initial part of the rod that is 
kept unmolten. Mostly, the zone passage is carried out upwards, because in that case the process has a higher stability. 

IS Usually, the growing crystal and sometimes also the melting rod are rotated, mostly with different rotation rates or coun- 
terrotation. 

The shape and stability of the molten zone play an important role in floating zone melting. The melt is held against 
gravity essentially by surface tension. In general, a molten zone Is more stable the higher its surface tension Is and the 
lower its specific weight. Besides, additional forces can act on the melt. In particular, when a high frequency induction 
20 heating system is used In semiconductor growth, skin effect generates non-negligible surface forces which help holding 
the molten zone up, while increasing melt convection. 

Aqueous solution, flux, and hydrothermal growth: 

25 The aqueous solution growth technique can be used to grow several salt single crystals. A simple aqueous solution 
crystallizer commonly consists of a rotating shaft with perpendicular branches, with a seed crystal attached to each 
arm. The fluid flow so provided breaks up the boundary layer of rejected solvent, leads to somewhat faster growth, and 
allows one to obtain more perfect crystals than those grown In quiescent solutions. Supersaturation is maintained by 
slow cooling, the rates being indeed very slow (0.1 to 1 °C per day). This requires a very careful and reliable tempera- 

30 ture control. 

The flux growth process is analogous to crystal growth from aqueous solutions, but the solvent solidifies before 
reaching room temperature. The main advantage of this method is that crystals are grown below the melting tempera- 
ture. If the material melts incongruently. I.e. decomposes before melting, or exhibits a phase transition below the melting 
point, one has indeed to look for growth temperatures lower than these phase transitions. Flux growth might thus be 

35 used when the melting temperature Is very high, and is useful when the vapor pressure at the melting temperature is 
too high. Thermal strain is minimized, due to the relatively low growth temperature, the very small temperature gradi- 
ents, and the free growth into a liquid, allowing the formation of growth facets. The method is very versatile, since a sol- 
vent may be found for any material required in the form of a single crystal. Flux growth is also suitable for the growth of 
layers on single crystal substrates. The main disadvantages are the low growth rate, the faceted crystal form and the 

40 unavoidable presence of ions of the flux as impurities in the crystals, if the solvent contains additional elements. In gen- 
eral, crystal growth from the melt Is preferable whenever possible. However, since often only small crystals are needed 
for basic investigations, the effort might be much smaller than that of melt techniques. 

Some very large crystals can also be grown by hydrothermal synthesis. This technique consists in using aqueous 
solvents or mineralizers under high temperature and pressure, in order to dissolve and recrystallize materials that are 

45 relatively Insoluble under ordinary conditions. Usually, In hydrothermal growth, dissolution Is carried out at the lower hot- 
ter zone (nutrient zone) and crystallization is carried out at the upper cooler zone (growth zone). Natural convection cur- 
rents created by the temperature gradient carry the material from the nutrient zone to the growth zone. Note that there 
are hydrothermal experiments carried out under extreme pressure conditions (less than 1 bar and greater than 10 
kbar). 

50 

Physical vapor transport growth: 

Bulk crystal growth can sometimes be performed by physical vapor transport. The advantage of these techniques 
Is that crystals tend to have a low concentration of point defects, and low dislocation densities, compared to crystals 
55 grown from the melt, the reason being that temperatures are usually considerably lower than the melting temperature. 
Moreover, If the material undergoes a phase transformation, or if it melts incongruently, vapor growth can be the only 
choice for the growth of single crystals. Compared to melt growth, a disadvantage of vapor growth Is however the rela- 
tively low growth rate typically encountered. 
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Several classes of techniques exist for the growth of bulk single crystals from the vapor phase. They are differenti- 
ated by the nature of the source material (e.g.. whether it is solid or vapor) and by the mechanisms by which it is trans- 
ported to the growing crystal surface. The conceptually simplest technique is that of sublimation, where the source 
material is placed at one end of a sealed tube and heated, so that It sublimes and is then transported to the colder 
5 region of the tube where it crystallizes. If the source material is a dissociable compound, then we have dissociative sub- 
limation, the transport processes are more complex and crystal stoichiometry becomes an issue to be considered. 

To further control the transport processes, an inert carrier gas can be introduced into the sealed tube. Transport 
can be facilitated by means of a reversible chemical reaction, where an active vapor species (often a halogen) reacts 
with the charge to produce a volatile species (a halide), which is transported to the crystallization zone. There, the reac- 
ts tion is reversed, leading to the deposition of the element or compound to be crystallized with release of the carrier (hal- 
ogen), which is recycled. 

Epitaxy: 

75 Thin film deposition processes are related to the bulk growth techniques in many aspects. In summary, one or sev- 
eral very thin layers of material are deposited onto a crystalline substrate to form a layered crystalline film. In epitaxial 
growth, many methods deal with Continuum Mechanics, namely Liquid and Solid Phase Epitaxy (including Halogen 
Transport Epitaxy, Levitation Epitaxy and Organo-Metallic Vapor Phase Epitaxy). All these techniques are very slow and 
concern very high quality products. 

20 

Factors affecting crystal quality: 

Regardless which growth technique is considered, the crystal quality is affected by three factors, namely the regu- 
larity of the lattice, the presence of impurities in the product, and the inhomogeneity of the crystal stoichiometry. 

25 While lattice defects can take the form of interstitials and vacancies, which are in part generated along the solidifi- 
cation front, and diffuse and recombine according to a complex thermally dependent kinetics, dislocations (or other 
irregularities, such as twins, grain boundaries etc.) are often induced by the thermal stresses which take place in the 
crystal near the solid-liquid interface. A key ingredient for obtaining high quality crystals is therefore to accurately control 
the evolution of the temperature field in the solid during and after growth. This can be achieved by regulating the thermal 

30 environment of the growing crystal and. in particular, by using radiative heat shields. 

On the other hand, impurities can be incorporated into the crystal during the growth process, either after having 
been transported through the melt from the surrounding atmosphere and container (crucible or ampoule), or as a con- 
sequence of a non perfect initial composition of the raw material. 

The third problem, which can occur when compounds or alloys are grown, is segregation. When several ions or 

35 molecules are present in the melt, the composition of the grown crystal can vary axially and radially, due to a preferential 
incorporation of one constituent in the solid phase. Segregation is mostly governed by the melt flow and, to a much 
lesser extent, by thermal effects, since the melting temperature can depend on crystal and melt composition. Obtaining, 
a uniform crystal composition is of course a key objective when e.g. lll-V compounds (such as GaAs crystals) are 
grown. 

40 Although any growth process involves peculiar difficulties, these are in general quite similar in nature and can be related 
to the three main quality factors hitherto described. 

The key issue considering the Czochralski process (Cz process) is to produce crystals whose diameter is almost 
constant and whose quality is the highest possible. This is a difficult objective, since most relevant physical quantities 
(such as interface deflection, growth speed of the lattice etc.). which characterize crystal quality, cannot at all be meas- 

45 ured during growth, in view of 

the high temperatures prevailing during the process, 

the extreme sensitivity of the crystal lattice with respect to any change of the environment of the liquid-solid inter- 
face. 

50 • the inherent batch nature of the process, and 
the high sensitivity to any contamination. 

Silicon growth: 

55 In particular, the quality of silicon crystals is specified according to the requirements of the customers. Such spec- 
ifications for the bulk quality include the concentration of oxygen, the resistivity, the density of micro defects and the 
maximum concentrations of metallic impurities, e.g. Fe. Cu. A case of particular interest is the problem of oxygen incor- 
poration in silicon crystals grown by the Czochralski process. As only quartz crucibles can be used, the solid SiOa is 
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dissolved partly into the melt and is evaporated mostly from the melt-gas interface. A fraction of the oxygen remains 
however dissolved in the liquid and arrives at the crystal-melt Interface, where a significant amount is Incorporated into 
the solid. The presence of oxygen in silicon crystals gives rise to positive and negative quality aspects. Oxygen precip- 
itates (8MD) in the center of the wafer getter metallic impurities during device manufacturing, which reduces current 
5 leakage. The dissolved oxygen in the wafer hardens the silicon lattice and makes it less susceptible for geometrical 
deformations. On the other side, BMD are detrimental for the devices, if they are located near the wafer surface. In this 
respect, the objective is to obtain a very homogeneous and reproducible oxygen content in the crystal within tight spec- 
ifications. 

The concentration of the dopants, such as B, P, As, Sb, determines the resistivity of the wafers. In order to obtain a 
10 high yield in device manufacturing, tightly specified ranges of dopant concentrations have to be incorporated homoge- 
neously into the crystal. 

However, during growth, the concentration of oxygen and dopants are non-measurable parameters, which signifi- 
cantly determine the quality of silicon crystals. Usually, the optimization of oxygen and dopant incorporation into the 
crystal is done by trial and error. Therefore, the development of crystal growth processes is very time consuming and 
15 expensive. 

During growth, so-called as-grown defects form, including vacancies, self-lnterstltials. oxygen and other impurity 
atoms, which diffuse and react with each other according to complex thermally dependent kinetics. A key ingredient for 
obtaining high quality crystals isthus again to accurately control the evolution of the temperature field in the solid during 
and after growth. This can be achieved by regulating the thermal environment of the growing crystal. In Czochralski sll- 
20 icon growth, peculiar components act as radiative heat shields, while the pull rate is selected in such a way that the ther- 
mal history of the crystal fits with the best experimental results. 

Measurement techniqes: 

25 Temperature measurements are very difficult to perform at high temperatures, since thermocouples can easily 
experience fast degradation, loss of accuracy etc., while pyrometers are In principle less accurate and more difficult to 
calibrate. Any on-line intrusive measurement in the melt (or the gas. In physical vapor growth) is strongly prohibited, 
since melt (or gas) composition cannot at all be modified during growth, in view of the severe composition requirements 
of most single crystals In high-tech applications. 

30 

Possible control strategies: 

A single control strategy is all the more not feasible, since the system geometry can undergo quite important 
changes during a full growth. In particular, the Gz process can be subdivided in several stages: melt down and stabili- 
35 zation, seeding, cone growth, roll-over or shouldering, body growth, tail-end growth and after-heating. An optimal con- 
trot structure would therefore require the use of a system of controllers that can continually adapt to the evolution of the 
geometry, as dictated by the crystal length increase and the melt volume decrease. Besides these difficulties, the 
growth control strategy is however facilitated in most cases by the process slowness. (Growing semi-conductor crystals 
can take several hours, or even a full day. Growing oxide crystals can take several days, or even weeks.) 

40 

Off-line control: 

In Czochralski growth, present control structures are able to control the diameter or the weight of the solidified crys- 
tal, which are the only easily measurable quantities, by acting on the pull rate and the heater power. No on-line quality 

45 control is presently available, and quality can only be controlled by off-line numerical simulations (with a view to adapt- 
ing the processing conditions). Finite element codes, such as 'FEMAG', are devoted to perform off-line simulations 
(FDupret, PNicodeme, Y Ryckmans, RWouters. M.J.Crochet, IntJ.Heat Mass Transfer. 33 (1990) 1849; M.J. Crochet. 
F Dupret, Y. Ryckmans. FT. Geyling; E.M, Monberg, J.Cryst. Growth, 97 (1989) i 73; F Dupret and N.van den Bogaert. 
in: Handbook of Crystal Growth. Vol. 2b, Chapter 15, edited by D.TJ. Hurle, North, Holland, 1994; R. Assaker, N. van 

50 den Bogaert. F Dupret. "Dynamic global simulation of bulk crystal growth under the influence of an axisymmetric mag- 
netic field", proceedings of the second international conference on energy transfer in magnetohydrodynamic flows, Aus- 
sois, France, 1994) which are employed to predict the temperature distribution and the shape of the free boundaries in 
all parts of a crystal puller, including crystal, melt, crucible, graphite susceptor. heater and several insulators as a func- 
tion of time, if necessary. The reference discusses in detail the underlying physical and numerical principles of the finite 

55 element code 'FEMAG'. which are summarized In this section. Based on the assumption of an axial symmetric furnace 
geometry, the code takes into account heat transfer by conduction in solid parts (e.g. crystal, heater, insulation) and 
radiation in the furnace enclosures. Convective heat transfer in the melt is computed by solving the complete Navier- 
Stokes equations with appropriate turbulence models. To reduce the computation time. It is possible as well to approx- 
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imate the convective heat transfer in the melt by means of an enhanced equivalent thermal conductivity. Heat transport 
by inert gas convection is usually neglected in silicon growth due to low pressure (ca. 20 mbar). This effect can be mod- 
elled approximately, when necessary. 

The shape of the solid/liquid interface, which is a priori unknown, is determined by the heat flux balance In the vicin- 

5 ity of the interface. That is. the heat flux q^ from the melt to the interface and the heat generated by solidification qf is 
equal to the heat flux qc into the crystal. Additionally the temperature at the solid/liquid interface Tgi is imposed to be 
equal to the melting temperature T^. 

In the furnace enclosures, heat is transferred by radiation, which couples the heat fluxes that are emitted, absorbed 
and reflected by solid and liquid surfaces. While all surfaces are usually assumed to be gray emitters, a band energy 

10 approach can be used to model wavelength dependent properties. This is necessary for semi-transparent materials 
(such as encapsulant layers or quartz tubes), which are modelled as transparent in a given range of frequencies and 
as opaque in the remaining of the spectrum. Diffuse radiation is taken into account, while specular radiation is 
neglected. A modified Gebhard method was applied to calculate the incoming and outgoing heat fluxes from the sur- 
faces. View factors are computed by taking viewed and hidden surfaces into account. 

15 The heat transfer equations form a set of non-linear, coupled, axisymmetric two-dimensional partial-differential 
equations, which are solved numerically by the finite element method. Galerkin's method is applied for space discreti- 
zation and the resulting set of nonlinear algebraic equations is solved by Newton's method. 

Data of the physical properties of all materials used in the simulations are required and can be selected from liter- 
ature or from the suppliers of the consumables. Thermal boundary conditions at the furnace walls have to be defined. 

20 e.g. 320 K at the furnace walls. - 

Present control strategies in silicon growth: 

Each stage of the growth is generally controlled through a different strategy The body phase regulation can include 
25 a system with two controllers, which regulate the diameter, on the one hand (by acting on the pull rate through a PD 
controller), and the average pull rate, on the other hand (by acting on the heater power through a PI controller). This 
regulation strategy for body growth is most often satisfactory with regard to the required set points for diameter or pull 
rate. However, control of the body phase for very large crystals is likely to represent a new difficulty, since, in that case, 
the process during body phase cannot be considered anymore as stationary and might become unstable due to large 
30 time constants and response times of the power controller. 

The difficulties of existing control strategies mainly deal with the regulation of the other stages of the process. Dur- 
ing seeding and cone growth, partial control only is achieved. In particular, a predetermined heater power evolution can 
be applied, while diameter is controlled via the pull rate. Control of the stages of melt down, roll-over and tail-end is usu- 
ally achieved by applying off-line precalculated set-point curves for heater power and pull rate. The control of the tail- 
35 end stage is a particularly difficult problem, since diameter cannot be measured today at that stage of growth. 

Other drawbacks of the current control strategies are that crystal quality is not directly addressed on-line, but only 
indirectly through actions on the pull rate or the diameter 

Summary Qf The Invention 

40 

A much better result could be achieved through an on-line control of the real quality criteria (impurities, stresses, 
etc.). As the growth process is usually very slow for most materials, techniques based on on-line calculations, which will 
provide the evolution of critical non-measurable physical quantities as growth proceeds, can be considered. 

An objective of the invention is to develop an appropriate control system in order to improve crystal quality in single 
45 crystal growth processes and to meet the more and more severe requirements that apply to crystals to be used in elec- 
tronics, optics, etc. 

The optimization of crystal quality is to be achieved by a direct control of the most relevant physical quantities which 
govern this quality. Functional dependencies of quality parameters on growth parameters are, e.g. in Czochralski silicon 
growth macroscopic and microscopic axial and radial variations of the oxygen and dopant concentrations as functions 
50 of pull rate, crystal rotation, crucible rotation, inert gas flow, melt level, heating system, etc. 

However, since these physical quantities can. most often, not be measured during growth, models which allow for the 
evaluation of such variables in real time have been developed. 

Another objective of the invention is to use multivariate predictive control strategies, based on the use of on-line 
simulations of the above mentioned models in the control loop, in order to be able to achieve accurate control of process 
55 variables (including non-measurable variables). 

The invention provides a process for controlling the growth of a crystal, comprising: 

establishing an on-line simulation software working with a reduced number of variables, the reduction of variables 
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being performed by using a projection algorithm; 

speeding up the on-line simulation software by generating a data bank in which values of off-line precalculated var- 
iables are stored; 

5 

tuning the on-line simulation software by adjusting the results predicted by on-line simulations to the results 
obtained by off-line simulations and by measurements; and 

establishing a control loop and controlling at least one of the variables in real time, the control loop using the 
10 speeded up and tuned on-line simulation software as an on-line observer. 

The invention also provides an apparatus for controlling the growth of a crystal, comprising: 

a) a crystal pulling means for pulling a single crystal from a melt; 
75 b) process control means for controlling the crystal pulling means; 

c) on-line simulation means for providing the process control means with on-line predicted values of parameters 
influencing the pulling of the crystal; 

d) data bank means for speeding up the on-line simulation means by providing data of on-line predicted values and 
off-line calculated values of parameters influencing the pulling of the crystal; and 

20 e) off-line simulation mean? for providing the control means, the on-line simulation means and the data bank means 
with off-line calculated values of parameters influencing the pulling of the crystal. 

Brief Description Of The Drawings 

25 Fig.1 refers to the on-line tuning of a reduced model which is based on a detailed distributed parameters model. 
Fig.2 refers to an apparatus which is suitable to perform the inventive process. 

Detailed Descriotion Of The Invention 

30 The on-line simulations are based on a distributed parameters model (a niodel based on partial differential equa- 
tions describing the major physical laws governing the process). These simulations can predict the evolution of some 
of the process variables and these predictions will be used in a predictive control scheme for better control of the proc- 
ess. 

The main ideas of this invention can be applied to any bulk crystal growth process, such as Czochralski. liquid 
35 encapsulated Czochralski. Bridgman, float zone or even any kind of epitaxial growth processes. The main common fea- 
tures of these processes are; 

The processes are slow, with characteristic time constants of the order of 10 minutes to several hours. 
The models for these processes are usually based on several PDEs describing the heat and mass exchanges and 
40 the interface shapes in various regimes of thiB apparatus. 

• The models are complex, due to the nature of the problem itself (crystallization of a liquid phase), the various com- 
ponents that need to be modelled, the influence of convection, the presence of boundaries, etc. 
The main variables which need to be controlled accurately in order to improve the quality of the crystal (and even 
the stability of the process itself when large crystals are grown) are non -measurable. 

45 

A successful implementation of this invention relies on four main steps, which are here described in the case of the 
Czochralski growth of silicon, but are also valid for other materials and techniques: 

A reduced model (a simplified version) of a detailed off-line simulation software (for Czochralski growth, this is. for 
50 example, the FEMAG software) is established, in order to speed-up the computations by reducing significantly the 
number of unknowns characterizing the thermal environment and the geometry of the crystal and the melt, so that 
this reduced model can be used on-line. This means that simulations run much faster than the real process. This 
reduction is performed using a projection algorithm. Some parameters governing the reduced model are tuned by 
fitting the results with detailed simulations. 
55 • In order to increase the computation speed even more, data banks generated by off-line simulations using the 
detailed simulation software are used to store values for the radiative matrices and for the melt flow stream function, 
for different geometries. During on-line simulation, interpolation between selected radiative matrices and stream 
functions are performed as a function of the current growth parameters. 
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For the reduced model to be used on-line as an observer for non-measurable variables, a tuning is necessary. This 
tuning Is done both off-line, by comparison between the results of the detailed and reduced models, and on-line, in 
real time, during the actual growth of the crystal, by comparison between computed and measured values of some 
of the process variables. Some material constants will then be adjusted on-line so that the simulated and predicted 

5 variables be close to their correct (i.e. measured) values. 

Once the reduced model operates on-line as a reliable observer of the process variables, it is used to generate pre- 
dictions of those variables of the process which need to be controlled. This can include measured variables, as for 
instance the pull rate which is necessary in order to obtain a specified diameter, but also unmeasured variables 
which are strongly linked to quality of the crystal, as for example interface shape or temperature field on the wall of 

10 the crucible. Multivariable predictive control techniques, based on the minimization of a performance index taking 
into account the difference between predicted and desired values for the controlled process variables, will then be 
applied to control the real process. 

These four steps are described in more detail below. The description is based on the Czochralski growth process. 

15 For the complex system of crystal growth, different simulation models (e,g. FEMAG) have been developed using 
finite element or finite volume techniques. These models are detailed global models, since the various heat transfer 
modes are taken into account, including conduction and radiation within and between the constituents, convection in 
the melt and release of solidification heat at the liquid-solid interface. The equations describing the process are nonlin- 
ear partial differential equations or integral equations. In order to solve these equations, one relies on a fine space dis- 

20 cretisation (commonly using finite elements) together with time discretisation (e.g. using an implicit Euler method). The 
reduced model, aimed at reducing significantly the computing time, is based on a projection technique. 
An easy description of the system is obtained by subdividing the furnace into a set of macro-elements, which can be 
"radiative" (the furnace enclosures), "two-dimensional" (e.g.. the solid bodies), or "one-dimensional" (the thin shells). 
The melt-crystal is a particular two-dimensional macro-element, whose shape is a priori unknown, and where latent 

25 heat is released at the solid-liquid interface. The set of interfaces connecting the macro-elements constitutes the skel- 
eton of the global domain. At any new time step, each macro-element includes a set of unknowns, which are the new 
nodal temperatures, represented by the vector T. the incremental generalized nodal heat fluxes on the skeleton, repre- 
sented by the vector Q"^, and the geometrical unknowns, represented by nodal positions X. The nodal temperatures can 
be split into two parts, namely T^. which represents the temperatures on the skeleton, and T"^, which represents all 
- 30 other temperatures. 

The objective of the reduced model is to build up a system of equations whose number of unknowns is much lower 
than that of the detailed model. For this purpose, the detailed discrete system is projected onto a reduced solution sub- 
space: 



40 in which U"^ is the vector of reduced nodal temperatures on the skeleton. The size n of U"^ is much smaller than the size 
of T* (n « N). The matrix is rectangular and depends on X (describing the positions of nodal unknowns). Reduc- 
tions are defined considering each side of the skeleton separately. Typical simplifications consist in assuming a con- 
stant, piecewise linear, or quadratic temperature profile along a given side. 

The global vector of unknowns corresponding to the detailed model. Z = [T*,T*.X.W]^ , is thus reduced into 

45 U = [U "*".T*,X,W] ^ . The relationship between U and Z can be written as 

Z = DU, 

where D is defined as follows: 

50 

D jj = D* jj if i ^ N and j ^ n. 
= 5jj i > N or j > n . 

55 Note that internal nodal temperatures In the furnace and geometrical unknowns are not reduced, since the idea is 
to maintain all details for what concerns the critical region (melt and crystal), while allowing simplifications in the other 
parts of the furnace. In order to keep the balance between equations and unknowns, the reduction in unknowns is per- 
formed together with a projection of the equations, which is basically achieved by multiplying the equations by the matrix 
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D^. The detailed discrete system of equations, which can be written in the following compact form at any time step: 

A(X)Z + B(X)F(X)T'* = E(X) . 

5 is, therefore, reduced in the following way : 

d'^A(X)DU + d'^B(X)F(X)(D*U'')'^ = D^E(X) . 

Hence, the system size is drastically decreased. There is no more need to use a decoupled scheme for finding the 
10 solution, and a coupled Newton-Raphson method is used. 

The factor (D'*'U*)^, which needs to be calculated for every iteration since the vector U"^ is modified, induces an 
important computing time in the Jacobian evaluation (through matrix multiplications involving as first factor a submatrix 
of size (n x N) of d''"B(X)F(X), and as second factor the matrix 4{D"^U'*")^D*, of size (N x n). A simple solution to decrease 
this cost is to define the following approximation : 



75 



35 



which can be shown to be satisfactory as long as the reduced nodes are not too sparsely arranged on the skeleton. 
With this approximation, it becomes necessary to only perform the product of a submatrix of size (n x n) of 
20 D^B(X)F(X)D'^ by the vector 4(4J'^)^, of size n, which is much cheaper. 
Further improvements lie in the following points : 

instead of using the previous approximation, the following alternative approach can be introduced : 

25 (D^U*)"* = (D D t;*)^ = D{b uY 

in which a different vector of reduced unknowns. D = D U, \s used for the radiative sides; the size of this vector is n^, 
with n < ng « N; this approximation is more accurate than before; 

30 • in order to ensure that essential boundary conditions, which are imposed along the external wall of the furnace, be 
exactly (and not approximately) verified on the reduced nodes, the projection matrix can be slightly modified into 
another matrix D* respecting these conditions. 



With such improvements, the system can be written in the form: 

D*A{X)DU+ D*B{X)F{X)D{b U*)^ = D*E{X). 



The reduced model described so far can be shown to be very accurate, provided node arrangement and reduced 
interpolation functions be selected in order that reduged and detailed simulations fit well together. 

40 Besides reducing the size of the system through a projection onto a reduced solution subspace. an important ingre- 
dient of our method is to avoid expensive calculations during on-line simulations by storing some information, which is 
computed beforehand, into a data bank, and access it as necessary. Two main classes of information are stored in this 
way, namely the radiative matrices BF and the stream function characterizing the melt flow. Radiative matrices are pre- 
calculated for different geometries, i.e. for different combinations of crystal length, and previous and current growth 

45 angles. For every crystal length, a small number of values for the current and previous growth angles are considered. 
Moreover, the radial positions of the nodes defining the crystal wall are imposed to be the prescribed crystal radius, 
except for the trijunction node and a few nodes located just above it, whose radial positions remain free. The error intro- 
duced by such an approximation Is negligible, for the crystal radius is to be maintained constant in the process, and the 
detailed shape would anyhow be very close to a cylinder. Storing pre-calculated radiative matrices allows to save com- 

50 puter time; moreover, the method is also cheap in terms of disk space, since the information which is stored is the pro- 
jected matrix, (D"^)'''BFD+, of size (n x n) (or (D^fSFV, of size (n x ng)), which is much less than the original matrix BF. 
A similar method is applied for the melt flow: the Stokes stream function is pre-computed by means of quasi-steady sim- 
ulations for different combinations of melt height, crystal and crucible rotation. rates, pull rate or any other parameter 
(such as magnetic field strength), and this information is stored in the data bank. During real-time simulation, interpo- 

55 lation between selected radiative matrices and stream functions are performed as a function of the current growth 
parameters. 

Once the reduced model has been designed with computation times compatible with' real time control of the proc- 
ess, there are different ways to introduce it inside a control loop. The possible options share a common feature. One 
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needs to be sure that the variables which are predicted or estimated by the reduced model do not drift away from their 
actual values on the process itself, although these actual values are usually not measurable. An online tuning mecha- 
nism is hence necessary, which makes sure that the reduced model behaves similarly as the actual process when the 
same input variables (pulling speed, heating power, etc.) are applied to both. Since the equations describing the proc- 

5 ess and which are integrated in the reduced model depend on several parameters (essentially material constants, 
which are not always accurately known), the adjustment of the reduced model is performed on some of these parame- 
ters, and based on deviations between computed and measured values for some process output variables, as, e.g.. the 
diameter and one or two temperature measurements. Rgure 1 illustrates the online tuning mechanism which trans- 
forms the reduced model into an observer for internal process variables. The reduced model is running on-line, using 

10 actual values for some of the input variables. For instance, it can use the actual power and pull rate, measured on the 
real process, to compute the solution of the PDEs. This solution will give the behavior of all other process variables. 
Some of these process variables can be measured on the real process and their measured values compared to the 
computed ones. Some parameters of the reduced model, typically some poorly known material constants, are then 
adjusted on-line in such a way that the difference between computed and measured values decreases with time. Differ- 

75 ent observer structures can be designed, according to which variables are chosen to act as inputs of the reduced 
model, which variables are selected as tuning variables and which parameters are used for the tuning mechanism. 

In all cases, once the reduced model works well, so that the difference between measured and computed variables 
is small enough, the other internal variables, which are computed by the reduced model but cannot be measured, can 
be used in a control loop. 

20 The internal variables whose future behavior can be predicted by the reduced model, are used in a multivariable 
predictive control loop. Crystal quality requirements are first transformed into set-point curves for some of these process 
variables. Then, a controller is designed which minimizes, over the set of admissible trajectories for the control varia- 
bles, the mean squared future deviations between the set-point curves and the predicted values. The criterion to be 
minimized may also include a second term, which takes into account the control energy which is required to achieve the 

25 objective. Typically, the criterion can be expressed in the following form: 

J(t) = Y^Mt+i) - y*(f+/)) ''oCyU+n ' + S Au(f+/) ''RAu{t+i) 

30 

where y (t+i) is a vector containing the predicted values of the process variables to be controlled at time t+i.y*(t+i) is 
the vector of set values for these variables and Au(t+i) is the vector of control increments. The prediction horizon Ny, the 
control horizon, N^, and the weighting matrices Q and R are tuning parameters for the general control algorithm. 
35 Since the model is highly nonlinear, minimization can be performed in at least three different ways: 

A nonlinear optimization algorithm can be used to search for a global minimum of the criterion. This way may be 
excessively costly in terms of computing time. 

One can compare the values of the criterion for different possible control trajectories and select the best one. 
40 • A satisfying control trajectory can be constructed on line by coupling the reduced model, working as a predictor, 
with a controller in a closed loop simulation. - v 

Depending on the global quality control objectives, several different internal variables may be selected as variables 
to be controlled and as control variables. 

45 With regard to the silicon Czochralski growth technique, measurable and non-measurable variables are preferably 
selected from a group comprising pulling speed, heater power, temperature distribution in the crystal, melt flow pattern, 
temperature distribution at the inner surface of the quartz crucible, temperature distribution on the melt surface, 
solid/liquid interface shape and evaporation of SiO. The oxygen and dopant concentration are controlled in such a way 
that the microscopic and macroscopic concentration variations in axial and radial direction of the crystal are small. 

50 Moreover, the axial oxygen content is preferably influenced by on-line controlling (some) non-measurable parameters, 
such as the temperature distribution at the inner surface of the quartz crucible (Si02 dissolution) and the evaporation 
rate of SiO from the melt surface. Finally, it is preferred to control the cooling of the growing crystal and in order to limit 
as-grown defects and to satisfy required specifications. 

Fig.2 refers to an apparatus which is suitable to perform the process for controlling the growth of a crystal. The 

55 apparatus comprises five units, which are connected. An off-line simulation unit predicts optimized process parameters, 
such as heater power, pull rate and rotation speeds of crystal and crucible (if, as an example, a Czochralski process is 
considered). These parameters are introduced into the process control unit. The on-line simulation control unit receives 
the same parameter set from the off-line simulation. 



10 
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The crystal growth process (crystal puller unit) Is controlled by the process control unit. Measurable parameters are 
fed back from the puller unit to the process control unit. The on-line simulation unit performs on-line simulations of the 
growth process in parallel. The on-line simulation unit receives the measurable parameters from the process control unit 
and compares them with the predicted values. The on-line simulation control unit sends new, predicted values for the 
5 future growth to the process control unit. The process control unit controls the crystal puller unit with the new values. 

A data bank unit is required to speed up the on-line simulation control unit. Process data from the off-line simulation 
unit and from the on-line simulation control unit are stored in the data bank where they can be retrieved by the on-line 
simulation control unit. The data bank unit contains permanent data from the off-line simulation unit and data that are 
up-dated at certain times by the on-line simulation control unit. 

70 

Advantages Of The Invention 

On-line simulation based control, that include functional dependencies for non-measurable quality parameters with 
respect to growth parameters, can reduce considerably the cost of crystal growth development. Moreover, the run to run 
15 stability of a growth process can be improved by an on-line simulation based control. 

Claims 

1 . Process for controlling the growth of a crystal, which growth is governed by a set of measurable and non-measur- 
20 able variables, comprising- 

establishing an on-line simulation software working with a reduced number of variables, the reduction of vari- 
ables being performed by using a projection algorithm; 

25 - speeding up the on-line simulation software by generating data banks in which values of off-line precalculated 
" variables are stored; 

tuning the on-line simulation software by adjusting the results predicted by on-line simulations to the results 
obtained by off-line simulations and by measurements; and 

establishing a control loop and controlling at least one of the variables in real time, the control loop using the 
speeded up and tuned on-line simulation software as an on-line observer. 

Process as claimed in claim 1, wherein the measurable and non-measurable variables are selected from a group 
comprising pulling speed, heater power, temperature distribution in the crystal, melt flow pattern, temperature dis- 
tribution at the inner surface of the quartz crucible, temperature distribution on the melt surface, solid/liquid inter- 
face shape and evaporation of SiO. 

An apparatus of controlling the growth of a crystal, comprising: 

a) a crystal pulling means for pulling a single crystal from a melt; 

b) process control means for controlling the crystal pulling means; 

c) on-line simulation means for providing the process control means with on-line predicted values of parame- 
ters influencing the pulling of the crystal; 

d) data bank means for speeding up the on-line simulation means by providing data of on-line predicted values 
and off-line calculated values of parameters influencing the pulling of the crystal; and 

e) off-line simulation means for providing the control means, the on-line simulation means and the data bank 
means with off-line calculated values of parameters influencing the pulling of the crystal. 

50 
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